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ABSTRACT 


The interstellar medium is turbulent and this induces relative motions of dust 
grains. We calculate relative velocities of charged grains in a partially ionized 
magnetized gas. We account for anisotropy of magnetohydrodynamic (MHD) 
turbulence, grain coupling with magnetic field, and the turbulence cutoff arising 
from the ambipolar drag. We obtain grain velocities for turbulence with parame¬ 
ters consistent with those in HI and dark clouds. These velocities are smaller than 
those in earlier papers, where MHD effects were disregarded. Finally, we con¬ 
sider grain velocities arising from photoelectric emission, radiation pressure and 
H 2 thrust. These are still lower than relative velocities induced by turbulence. 
We conclude that turbulence should prevent these mechanisms from segregating 
grains by size. 

Subject headings: ISM:dust, extinction—kinematics,dynamics—magnetic fields 


1. Introduction 

Dust is an important constituent of the interstellar medium (ISM). It interferes with 
observations in the optical range, but provides an insight to star-formation activity through 
far-infrared radiation. It also enables molecular hydrogen formation and traces the magnetic 
field via emission and extinction polarization. The basic properties of dust (optical, align¬ 
ment etc.) strongly depend on its size distribution. The latter evolves as the result of grain 
collisions, whose frequency and consequences depend on grain relative velocities. 

Various processes can affect the velocities of dust grains. Radiation, ambipolar diffusion, 
and gravitational sedimentation all can bring about a dispersion in grain velocities. It 
was speculated in de Oliveira-Cost a et al. (2000) that starlight radiation can produce the 
segregation of different sized grains that is necessary to explain a poor correlation of the 
microwave and 100 pm signals of the foreground emission (Mukherjee et al. 2001). If true 
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it has big implications for the CMB foreground studies. However, the efficiency of this 
segregation depends on grain random velocities, which we study in this paper. 

Interstellar gas is turbulent (see Arons & Max 1975). Turbulence was invoked by a 
number of authors (see Kusaka et al. 1970, Volk et al. 1980, Draine 1985, Ossenkopf 1993, 
Weidenschilling & Ruzmaikina 1994) to provide substantial relative motions of dust particles. 
However, they discussed hydrodynamic turbulence. It is clear that this picture cannot be 
applicable to the magnetized ISM as the magnetic fields substantially affect fluid dynamics. 
Moreover dust grains are charged, and their interactions with magnetized turbulence is very 
different from the hydrodynamic case. This unsatisfactory situation motivates us to revisit 
the problem and calculate the grain relative motions in magnetized ISM. In what follows, we 
use the model of MHD turbulence by Goldreich and Sridhar (1995, henceforth GS95), which 
is supported by recent numerical simulations (Clio & Vishniac 2000, Maron & Goldreich 
2001, Clio, Lazarian & Vishniac 2002a, henceforth CLV02). We apply our results to the cold 
neutral medium (CNM) and a dark cloud to estimate the efficiency of coagulation, shattering 
and segregation of grains. 


2. MHD Turbulence and Grain Motion 

Unlike hydrodynamic turbulence, MHD turbulence is anisotropic, with eddies elongated 
along the magnetic held. This happens because it is easier to mix the magnetic held lines 
perpendicular to their direction rather than to bend them. The energy of eddies drops with 
the decrease of eddy size (e.g. Vi ~ l 1 ^ 3 for the Kolmogorov turbulence) and it becomes 
more difficult for smaller eddies to bend the magnetic held lines. Therefore the eddies get 
more and more anisotropic as their sizes decrease. As eddies mix the magnetic held lines 
at the rate k±Vk, where k±is a wavenumber measured in the direction perpendicular to the 
local magnetic held and ry is the mixing velocity at this scale, the magnetic perturbations 
propagate along the magnetic held lines at the rate k\\VA ,where k\\ is the parallel wavenumber 
and Va is the Alfven velocity. The corner stone of the GS95 model is a critical balance 
between those rates, i.e., k±Vk~ k\\VA, which may be also viewed as coupling of eddies and 
wave-like motions. Mixing motions perpendicular to the magnetic held lines are essentially 
hydrodynamic (see CLV02) and therefore it is not surprising that the GS95 predicted the 
Kolmogorov one-dimensional energy spectrum in terms of k±, i.e., E{k±) ~ A. - j_ 5 ^ 3 (see review 
by Clio, Lazarian & Yan 2002, henceforth CLY02). 

The GS95 model describes incompressible MHD turbulence. Recent research suggests 
that the scaling is approximately true for the dominant Alfvenic modes in a compressible 
medium with Mach numbers(M = VjC s ) of the order of unity (Lithwick & Goldreich 2001, 



3 


henceforth LG01, CLY02, Clio & Lazarian 2002, in preparation), which is also consistent 
with the analysis of observational data (Lazarian & Pogosyan 2000, Stanimirovic & Lazarian 
2001, CLY02). In what follows we apply the GS95 scaling to handle the problem of grain 
motions. 


Because of turbulence anisotropy, it is convenient to consider separately grain motions 
parallel and perpendicular to the magnetic field. The motions perpendicular to the magnetic 
field are influenced by Alfven modes, while those parallel to the magnetic field are subjected 

_1/3 

to the magnetosonic modes. The scaling relation for perpendicular motion is oc k ± 7 
(GS95). As the eddy turnover time is t*, oc (k^Vk) -1 , the velocity may be expressed as 
Vk ~ V max {Tk/T max ) 1/2 , where T max = lmax/v max is the time-scale for the largest eddies, for 
which we adopt the fiducial values l max = lOpc, v max = 5km/s. 


Grains are charged and coupled with the magnetic field. If the Larmor time tl = 
27T m gr c/qB is shorter than the gas drag time tdrag , grain perpendicular motions are con¬ 
strained by the magnetic field. In this case, grains have a velocity dispersion determined 
by the turbulence eddy whose turnover period is ~ tl, while grains move with the mag¬ 
netic field on longer time scales. Since the turbulence velocity grows with the eddy size, the 
largest velocity difference occurs on the largest scale where grains are still decoupled. Thus, 
following the approach in Draine (1985), we can estimate the characteristic grain velocity 
relative to the fluid as the velocity of the eddy with a turnover time equal to tl, 

" w (s?)'"■•" “> 

and the relative velocity of grains to each other should be approximately equal to the larger 
one of the grains’ velocities, i.e., the the larger grain’s velocity, 


5u_l(oi, a 2 ) 


7 i ^(^) i/2 (^) 1/2 [^(«.,« 2 )] 3/2 

lmax \ Q / 

1.4 x 10 5 cm/s(v 5 a 5 ) 3/2 / (q e l w B g ) 1/2 , 
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in which v 5 = n ma3 ,/10 5 crn/s, a 5 = a/10 5 cm, q e = g/lelectron, l 10 = l max / lOpc, B g = 
B/lpG, and the grain density is assumed to be p gr = 2.6g/cm -3 . 


Grain motions parallel to the magnetic field are induced by the compressive compo- 

_i /o 

nent of slow mode with v\\ oc k^ 7 (CLV02, LG01, CLY02). The eddy turnover time is 

Tfc oc (u||/c||) _1 , so the parallel velocity can be described as uy ~ v max Tk/T max x . For grain mo¬ 
tions parallel to the magnetic field the Larmor precession is unimportant and the gas-grain 


1 We assume that turbulence is driven isotropically at the scale l max - 
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coupling takes place on the translational drag time td rag • The drag time due to collisions 
with atoms is essentially the time for collision with the mass of gas equal to the mass of 
grain, t drag = (ap 9r /n)(7r/8/i/cT) 1//2 ., where // is the mass of gas species. The ion-grain cross- 
section due to long-range Coulomb force is larger than the atom-grain cross-section (Draine 
& Salpeter 1979). Therefore, in the presence of collisions with ions, the effective drag time 
decreases, tdrag = , where a < 1 is the function of a particular ISM phase. The char¬ 

acteristic velocity of grain motions along the magnetic held is approximately equal to the 
parallel turbulent velocity of eddies with turnover time equal to td rag 


U||(a) 


= a 


v 

max 

Imax 




and the relative velocity of grains for T 10 o = T/100K is 
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(1.0 x 10 6 cm/s)au 5 a 5 /(n/ 10 T 1 1 0 / 0 2 ), 
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When t j > tdrag, grains are no longer tied to the magnetic held. Since at a given scale, 
the largest velocity dispersion is perpendicular to the magnetic held direction, the velocity 
gradient over the grain mean free path is maximal in the direction perpendicular to the 
magnetic held direction. The corresponding scaling is analogous to the hydrodynamic case, 
which was discussed in Draine (1985): Sv(ai,a 2 ) = vUlx/lmax.t]/ r 2 ag . i.e., 


Sv(a i, a 2 ) 


3/2 
I Vmax 

a 2 - 


1/2 
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Turbulence is damped due to the viscousity when the cascading rate v±k± equals the 
damping time tdamp (see Clio, Lazarian & Vishniac 2002b). If the mean free path for a 
neutral particle l n , in a partially ionized gas with density n tot = ri n + n % , is much less 
than the size of the eddy in consideration, i.e., l n k± -C 1, the damping time is td amp ~ 
W'fcJ 2 ~ (n tot /n n ) (InVn )- 1 kj 2 , where u n is effective viscosity produced by neutrals. In the 
present paper we consider cold gas with low ionization, therefore the influence of ions on l n 
is disregarded. Thus the turbulence cutoff time in neutral medium is 
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where v n and Va are, respectively, the velocity of a neutral and Alfven velocity. It is easy to 
see that for r c longer than either td ra g or a the grain motions get modified. A grain samples 
only a part of the eddy before gaining the velocity of the ambient gas. In GS95 picture, 
the shear rate dvjdl increases with the decrease of eddy size. Thus for r c > max{td ra g,TL}, 
these smallest available eddies are the most important for grain acceleration. Consider first 
the perpendicular motions. If v c is the velocity of the critically damped eddy, the distance 
traveled by the grain is A/ ~ v c x mm{Grac/, a}- Thus the grain experiences the velocity 
difference A Ixdv/dl ~ v c xmin{td ra g , A,} Ac- Due to the critical balance in GS95 model, the 
shear rate along the magnetic field is dv/dl = v c k\\ ~ v c /(Vat c ). Therefore, grain experiences 
a velocity difference Va/v c times smaller, i.e., ~ v 2 c x td rag / (U 4 T C )- 


3. Discussion 

3.1. Shattering and Coagulation 

Consider the cold neutral medium (CNM) with temperature T = 100K, density n H = 
30cm -3 , electron density n e = 0.045cm -3 , magnetic held B ~ 1.3 x 10 -5 G (Weingartner 
& Draine 2001a, hereafter WDOla). To account for the Coulomb drag, we use the results 
by WDOla and get the modified drag time tdrag = a tdrag■ Using the electric potentials in 
Weingartner & Draine (2001b), we get grain charge and ry. 

For the parameters given above, we find that td ra g is larger than r c for grains larger than 
10 -6 cm, T/\ is smaller than r c even for grains as large as 10 -5 cm. Here, we only consider grains 
larger than 10 -6 cm, which carry most grain mass (~ 80%) in ISM, so we can still use Eq.(3) 
to calculate grain parallel velocities and Eq.(l) to get the perpendicular velocity for grain 
larger than 10 -5 cm. Nevertheless, the perpendicular velocities of grains smaller than 10 -5 cm 
should be estimated as t/^a) = v c x (t l /t c ) = v max (T c /T max ) 1/2 (T L /T c ) = v ± (o)(t l /t c ) 1/2 , 
where v±(a) is given by Eq.(l). The results are shown in Fig.l. 

The critical sticking velocity were calculated in Chokshi et al. (1993) (see also Dominik & 
Tielens 1997). 2 However, experimental work by Blum (2000) shows that the critical velocity 
is an order of magnitude larger than the theoretical calculation. Thus the collisions can 
result in coagulation for small silicate grains (< 3 x 10 -6 cm). 

With our input parameters, grains do not shatter if the shattering thresholds for silicate 


2 There are obvious misprints in the numerical coefficient of Eq.(7) in Chokshi et al.(1993) and the power 
index of Young’s modulus in Eq.(28) of Dominik & Tielens (1997). 
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is 2.7km/s as in Jones et al. (1996). Nevertheless, the grain velocities strongly depend on 
v max at the injection scale. For instance, we will get a cutoff 6 x 10 -5 cm due to shattering 
if v max = lOkm/s. 

For a dark cloud, the situation is different. As the density increases, the drag by 
gas becomes stronger. Consider a typical dark cloud with temperature T = 20K, density 
n h = 10 4 cm -3 (Chokshi et al. 1993) and magnetic held B ~ 2.3 x 10 _4 G. Assuming that 
dark clouds are shielded from radiation, grains get charged by collisions with electrons: 
< q >= 0.3(r/10 _5 cm) electrons. The ionization in the cloud is x — n e/ n tot ~ 10” 6 and the 
drag by neutral atoms is dominant. From Eq.(6) and the expression for the drag time and 
the Larmor time, we find tl < tdrag for grains of sizes between 10 _6 cm and 4 x 10 _6 cm, and 
tdrag < tl for grains larger than 4 x 10 _6 cm. In both cases, turbulence cutoff r c is smaller 
than tdrag and r L . Thus for the smaller grains, we use Eq.(l),(3) to estimate grain velocities. 
For larger grains, grain velocities are given by Eq.(5). 

Our results for dark clouds show only a slight difference from the earlier hydrodynamic 
estimates. Since the drag time td rag oc n -1 , Larmor time ry oc B~ l oc n -1 / 2 , the grain 
motions get less affected by the magnetic held as the cloud becomes denser. Thus we agree 
with Chokshi’s et al. (1993) conclusion that densities well in excess of 10 4 cm -3 are required 
for coagulation to occur. Shattering will not happen because the velocities are small, so 
there are more large grains in dark clouds. This agrees with observations (see Mathis 1990). 

In the treatment above we disregarded the possibility of direct acceleration of charged 
grains through their interactions with fluctuating magnetic held. In our next paper we will 
show that this resonant process is important for a highly ionized medium. 


3.2. Grain Segregation and Turbulent Mixing 

Our results are also relevant to grain segregation. Grains are the major carrier of heavy 
elements in the ISM. The issue of grain segregation may have signihcant inhuence on the ISM 
metallicity. Subjected to external forcing, e.g., due to radiation pressure, grains gain size- 
dependent velocities with respect to gas. WDOla have considered the forces on dust grains 
exposed to anisotropic interstellar radiation helds. They included photoelectric emission, 
photodesorption as well as radiation pressure, and calculated the drift velocity for grains of 
different sizes. The velocities they got for silicate grains in the CNM range from O.lcm/s to 
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10 3 cm/s. Fig.l shows that the turbulence produces larger velocity dispersions. 3 Thus the 
grain segregation of very small and large grains speculated in de Oliveira-Costa et ah (2000) 
is unlikely to happen for typical CNM conditions. 

A different mechanism of driving grain motions is a residual imbalance in “rocket thrust” 
between the opposite surfaces of a rotating grain (Purcell 1979). This mechanism can provide 
grain relative motions and preferentially move grains into molecular clouds. It is easy to see 
that due to averaging caused by grain rotation, the rocket thrust is parallel to the rotation 
axis. Three causes for the thrust were suggested by Purcell (1979): spatial variation of the 
accommodation coefficient for impinging atoms, photoelectric emission, and H 2 formation. 
The latter was shown to be the strongest among the three. The uncompensated force in this 
case arises from the difference of the number of catalytic active sites for H 2 formation on the 
opposite grain surfaces. The nascent H 2 molecules leave the active sites with kinetic energy 
E and the grain experiences a push in the opposite directions. The number of active sites 
varies from one grain to another, and we should deal with the expectation value of the force 
for a given distribution of active sites. 

Due to internal relaxation of energy (see Lazarian & Draine 1999a,b, and review by 
Lazarian 2000) the grain rotational axis tends to be perpendicular to the largest b — b 
surface. Adopting the approach in Lazarian & Draine (1997), we get the mean square root 
force of H 2 thrust on a grain in the shape of a square prism with dimensions bxbxa (b > a) 

(F zH ) = r 3/2 (r + 1) 1/2 7(1 - y)n n v n a 2 , ( 7 ) 

where r = b/2a„ rip = n(H) + 2n(H 2 ), y = 2n(H 2 )/nn is the H 2 fraction, 7 is the fraction 
of impinging H atoms and v is the number of active sites over the grain surface. The 
expected grain velocity is v = In the CNM we consider, y = 0, adopting the 

characteristic values in Lazarian & Draine (1997), r — 1,7 = 0.2, E = 0.2eV, and the density 
of active sites 10 11 cm -2 so that v = 80(a/10 _5 cm) 2 r(r + 1), we get the “optimistic” velocity 
shown in Fig 1. For maximal active site density 10 15 cm -2 , we get the lower boundary of 
grain velocity v ~ 3.3(10 _5 cm/a) 1 / 2 cm/s. The scaling is approximate due to the complexity 
of coefficient a(see WDOla Fig.16). 

Lazarian & Draine (1999a,b) have shown that subjected to H 2 torques alone, grains 
< 10 _4 cm should experience frequent thermal flipping, which means that the F z h fluctuates. 
This flipping results from coupling of grain rotational and vibrational degrees of freedom 


3 If reconnection is fast (see Lazarian & Vishniac 1999), the mixing of grains over large scales is provided 
by turbulent diffusivity~ v max l max . On small scales the grain decoupled motions are important. 



through internal relaxation and would average out (F~ h). However, the flipping rate depends 
on the value of the grain angular momentum (Lazarian & Draine 1999a). If a grain is already 
spun up to a sufficient velocity, it gets immune to thermal flipping. Radiative torques 
(Draine & Weingartner 1996) can provide efficient spin if the grain size is comparable to 
the wavelength. For a typical interstellar diffuse radiation field, the radiative torques are 
expected to spin up grains with sizes larger than ~4x 10 _6 cm. They will also align grains 
with rotational axes parallel to the magnetic field. Thus grains should acquire velocities 
along the magnetic field lines and the corresponding velocities should be compared with 
those arising from turbulent motions parallel to the magnetic field. It is clear from Fig.l 
that for the chosen set of parameters the effect of H 2 thrust is limited. All in all, we conclude 
that the radiation effects and H 2 thrust are not efficient for segregating grains in typical ISM 
conditions. 


4. Summary 

We have calculated relative motions of dust grains in a magnetized turbulent fluid 
taking into account turbulence anisotropy, turbulence damping and grain coupling with the 
magnetic field. We find that these effects decrease the relative velocities of dust grains 
compared to the earlier hydrodynamic-based calculations. The difference is substantial in 
CNM, but less important for dark clouds. For CNM we find that coagulations of silicate 
grains happen for sizes < 3 x 10 _6 cm. The force due to H 2 formation on grain surface 
might drive small grains (< 3 x 10 _6 cm) to larger velocities but thermal flipping of grains 
suppresses the forces for grains less than 4 x 10“ 6 cm. We conclude that radiation and H 2 
thrust are not capable of segregating grains. 

We are grateful to John Mathis for reading the manuscript and many important com¬ 
ments. We thank our referee Dr. Stuart Weidenschilling for helpful comments. The research 
is supported by the NSF grant AST0125544. 
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Fig. 1.— Grain velocities as a function of radii (solid line) in the CNM. Dashdot line 
represents parallel velocity due to the drag by compressible modes, dotted line refers to 
perpendicular velocity from the contribution of the drag by Alfven mode, also plotted is 
the earlier hydrodynamic result (dashed line). The change of the slope is due to the cutoff 
of turbulence by ambipolar diffusion. The grain velocity driven by H 2 thrust is plotted to 
illustrate the issue of grain segregation in the CNM (see text), the part marked by ’o’ is 
nonphysical because thermal flipping is not taken into account. 




